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Abstract 

Based on the spin-density-matrix approach, both the electric-field-induced spin accumulation and 
the spin current are systematically studied for the two-dimensional Rashba model. Eigenmodes of 
spin excitations give rise to resonances in the frequency domain. Utilizing a general and physically 
well-founded definition of the spin current, we obtain results that differ remarkably from previous 
findings. It is shown that there is a close relationship between the spin accumulation and the spin 
current, which is due to the prescription of a quasi-chemical potential and which does not result 
from a conservation law. Physical ambiguities are removed that plagued former approaches with 
respect to a spin-Hall current that is independent of the electric field. For the clean Rashba model, 
the intrinsic spin-Hall conductivity exhibits a logarithmic divergency in the low-frequency regime. 

PACS numbers: 72.25.-b 73.23.-b 73.50.Bk 



1 



I. INTRODUCTION 



It has been anticipated that the spin degree of freedom of charge carriers potentially 
provides additional functionality to electronic devices. Recent progress aims at spintronic 
applications that rely on the capability to manipulate electron spin polarizations in nonmag- 
netic semiconductors. Of particular interest are efficient injection mechanisms of spins in 
semiconductors at room temperature. In this respect, the proposal by Murakami et al. 
and Sinova et al. 2] of generating dissipationless transverse spin currents by a driving elec- 
tric field has attracted considerable attention. The observation of this spin-Hall effect has 
been reported in recent experiments on GaAs and related materials . Many theoretical 



studies of this effect 
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focused on the Rashba spin-orbit 



interaction since this type of spin-charge coupling can be easily controlled by an electrical 
gate. The original conclusion concerning the existence of an universal intrinsic spin-Hall 
conductivity ^ in clean systems was reexamined in more detail by treating effects of the 



elastic impuritv scattering. Based on numerical 



the Keldysh 



itv scati 

HQ 



and Kubo 



ITL 18^ and analytical results derived from 



2l| formalism, it has been concluded that vertex 



corrections lead to a vanishing zero-frequency spin-Hall current in the thermodynamic limit 
of the linear Rashba model. After long debates many researchers finally arrived at the same 
conclusion so that there seems to be agreement now that the intrinsic zero-frequency spin 
current is finite and universal for a free two-dimensional electron gas but vanishes in impure 
systems for an arbitrary ratio of spin splitting and the impurity scattering rate. 

Despite this consensus there are still challenging problems to be addressed referring to a 

)roper definition of the spin current. This issue has recently been treated by Zhang et al. 

22 1 and Sugimoto et al. 23|. The authors pointed out that the "conventional" spin current. 



which is defined as the product of spin and velocity operators, loses its physical foundation 
when the spin-orbit coupling is present. The main deficiency of this definition is related to 
the absence of a conservation law of spins. Therefore, a physically motivated definition was 
sugg ested HQ 

that relates the spin current to the time derivative of the spin displacement. 
It was argued that under quite general conditions this effective spin current satisfies the 
continuity equation so that it is measurable as a spin accumulation. While for the charge 
transport both definitions completely agree to each other, there is a remarkable discrepancy 
between them with respect to the spin current due to the torque dipole contribution j2z 
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Results for the spin Hall conductivity derived from this physically motivated definition 
of the spin current remarkably differ from previous findings based on the "conventional" 
definition. Unfortunately, the authors did not apply direct perturbational techniques to 
calculate this spin-Hall current. Rather, they worked out a special calculational schema 
for the determination of the torque dipole density, which they used to complement the 
"conventional" spin-Hall current to a conserved quantity. In view of the long and successful 



history to derive the current from the time derivative of the dipole moment [2J|, we prefer 
the application of the standard procedure. What is needed in this approach is nothing but 
the density matrix, which is determined from quantum-kinetic equations. Starting from 
the physical definition of the spin-Hall current, we calculate the frequency dependent spin 
polarization and spin-Hall conductivity by analytically solving the kinetic equations for the 
spin-density matrix. The obtained results do not agree with previous findings deduced from 
the "conventional" definition of the spin-Hall current. Furthermore, physical inconsistencies 
with respect to a spin-Hall current component that is independent of the electric field as 
well as the relationship between the spin accumulation and the spin current are puzzled out. 
For the spin-Hall conductivity of a clean two-dimensional electron gas with Rashba type 
spin-orbit coupling, our exact calculation does not reproduce any universal value but yields 
a logarithmic divergency in the low frequency limit. In addition, the frequency dependent 
spin-Hall conductivity, which is obtained by a controlled perturbational approach, differs 
even qualitatively from previous results [l^ . 

II. THE KINETIC EQUATIONS 

We consider a two-dimensional electron gas in the presence of Rashba spin-orbit inter- 
action with amplitude a. The system with an applied in-plane electric field E (which is 
oriented along the x axis) is described by the Hamiltonian 

fc,A fc,A,A' 

- eEj2^'^(^L^x(^k+^x\ , (1) 

' ' ^ I K.=0 



fc,A 

where we introduced the abbreviations 



£fc = — , = —{K xk), K = -TT^z- (2) 
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m, ef, and a denote the effective mass, the Fermi energy, and the vector of Pauh matrices, 
respectively, al.^ and akx are creation and annihilation operators with quasimomentum 
A; = {kj., ky, 0) and spin A. We are going to calculate the time-dependent density matrix 

f^,{k,k'\t) = {ai^ak'y)t, (3) 
which is more conveniently expressed by its physical representation 

f{k, i^\t) = J2 fxik, AC I t), f{k, i^\t) = J2 fx'ik, 1^ I t)^A,A', (4) 

A A,A' 

in the fc, n space, where k fc + K/2 and k' k — n/2. n refers to a possible inhomogeneity 
of the charge and/or spin distribution. Using the Liouville equation, the quantum-kinetic 
equations for the components of the density matrix are straightforwardly derived. From the 
result 

k)f - -K{f^ /.) + ^V./ (5) 

ot m m h 

= E E{/A'^(^>l^)<A\fc^fc,'^)-/A'^fc,'^l^)<\'(fe,fc^^)W^ 

A,Ai,A2 k' ^ ^ 



k)f + 2{Q X /) + -{K X + ^V./ (6) 

ot m m h 



\, \„ \., \, 1.1 ^ ' 



Ai,A2 A3,A4 k' 

it is concluded that not only scattering but also any inhomogeneity ^ 0) couples the 
charge (/) and spin (/) degrees of freedom to each other. Consequently, any accumulation 
of charges induces a spin response and vice versa. The left-hand side of these equations was 



derived and discussed in Ref. |25j. The scattering probabilities W^^^^ on the right hand 
sides of Eqs. (5) and (6) comprise both elastic and inelastic scattering-in and scattering- 
out contributions, which satisfy a sum rule. We shall restrict the consideration to elastic 
scattering described by the Hamiltonian 

Hint = " E E '^fcA'^fc'A; (7) 
k,k' A 

with u denoting the magnitude of the short-range impurity potential. For time- dependent 
phenomena, we prefer the treatment of the Laplace transformed kinetic Eqs. (5) and (6). 



Within the Bom approximation, we obtain for the scattering probabihties the exact result 



<\^^(fc',fc,K|s) = ^ / citexp 



— st + — (efc'-K/2 — £k+K./2) t 



X 



COS {uJk'-n/2t) 5aiA3 " 
COS {uJk+n/2t) 5xi\2 + i 



sm 



^k'-K./2 
<?A4A2 ■ '^fc+K/2 . 



(^fc'-fc/2^) 



sin (t^fc+^/2t) 



(8) 



^k+K,/2 

with s denoting the variable of the Laplace transformation. As usual, it is assumed that 
corrections due to the k expansion of the scattering probabilities are small compared to 
corresponding contributions on the left-hand side of the kinetic equations. Furthermore, for 
weak spin-orbit coupling, we may restrict to the lowest-order contributions in uj^t. Adopting 
these approximations, the collision integrals are expressed by 

1 d 



if - f) —i^^^kf + 



T 



T del ^ 



^(^k ■ f 



(9) 



(ek-ek'Y 



1 d 



hwu d 



f 



h 

+ - 

T 



del 



hwk X / 



2u 



fc' 



r dek 
OJk + l^k>) X /(fc') 



[^k — ^k' 



(10) 



UJk' X UJk 

{ey - SkY 



where the scattering time r is calculated from 

1 

k> 



h 



^5(£fc/ - Ek). 



(11) 



/(fc) means an average over the angle of the vector fc. In the Eqs. and ()10|) corrections 
appear, which result from virtual transitions that are not linked to the scattering time r. 
The treatment of these contributions as well as higher-order corrections in the spin-orbit 
coupling a goes beyond the scope of this paper. Restricting to the lowest-order scattering 
contributions, assuming an initial thermodynamic equilibrium state, and focusing on weak 
spin-orbit coupling so that h^K / {mek) <^ 1, the kinetic equations are expressed by 



sf - -{k ■ k)f - -Kifx A.) + ^Vfc/ = -(/ - /) + nisk) 
mm n T 



(12) 
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sf + 2(c3fc X /) - -{k ■ k)f + -{K X n)f + ^Vfc/ 

mm n 

= -{f - f) + -T^f^k —f-huk^ , (13) 

where n{ek) denotes the initial equihbrium charge density. 



III. GENERAL EXPRESSIONS FOR THE CHARGE AND SPIN CURRENTS 

The Laplace-transformed density matrix f{k,K, \ s) = {f{k,K, \ s),f{k,n \ s)} contains 
all the information needed to determine all kinetic observables. In particular, the quantity 

f{s) = J2fik,'^\s)U=o (14) 

k 

represents nothing but the conserved charge density. Exactly in the same way, the total 
magnetic moment is calculated from 

fc 

When an electric field is applied to the system, these quantities become time dependent and 
allow the treatment of relaxation processes. 

Other quantities of interest are the current of charge carriers and spins, which are obtained 
from 

= -ise^ V^/(fc, K I s) 1^=0, (16) 
fc 

and 

?(s) = -is^ ^ ® /(fc, K I s) 1^=0, (17) 
fc 

respectively. In Eq. (fTTj) . (g) denotes the dyadic product. These definitions are fundamental 
and sufficiently general. In the time domain, the spatial versions of these equations describe 
the temporal evolution of the carrier or spin displacement, i.e., the center-of-mass velocity of 
the wave packet. The expression for the charge- carrier current in Eq. is also applicable 
for systems without any spatial dispersion. If the interaction Hamiltonian that describes 
elastic or inelastic scattering commutes with the dipole operator of the carriers than the 
definition in Eq. (fTBj) becomes completely equivalent to 

^.w4Et^('=i')4E|i-(^x/('=U)),. (18) 

fc k,j ■' 
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in which the k dependence does no longer occur. We want to stress that the situation for 
the spin transport is completely different and more subtle. For a system without any spatial 
dispersion, it is in general not possible to express the spin current as defined in Eq. (|17p 
by the density matrix f{k,K, \ s) 1^=0 alone. What is really needed is the n derivative 
of this function at k = 0. Notwithstanding this fact, most researchers defined the spin 
current j^{s) as the symmetrized product of the spin and velocity operators {ai,Vk}+/4:, 
where v = VkHo/h and Ho is obtained from Eq. (P) for E = 0. For the Rashba model, this 
definition of the spin current takes the for 

" ^ 5Z ^'^^f' ® f^^' i^\s)\n=o. (19) 

k 

As shown below, the definition in Eq. (1191) leads to a stationary spin-Hall current that does 

nn 

not depend on the electric field. In addition, recent studies [22, |23|| clearly demonstrated 
that results derived from the Eqs. ()17|) and p9|) considerably differ from each other and 
that the measurable quantity is related to Eq. (fT7j) . Furthermore, we note that the diffusion 
tensor is likewise obtained from derivatives of the density matrix f{k,n \ s) with respect 
to at At = 0. Although the extension of our approach to the treatment of the diffusion 
coefficient is straightforward, we want to confine ourselves to the analysis of spin currents. 

The definitions of measurable quantities in Eqs. (14) to (18) already set up our calcula- 
tional schema as an iteration with respect to k. As we restrict ourselves to elastic scattering, 
the spectral functions can be calculated for a given energy over which one finally integrates. 
This procedure is applicable only in the linear response regime, where carrier heating and 
energy relaxation due to nonlinear field effects do not play an essential role. Therefore, we 
treat only first-order corrections in the electric field. The analytic solution of the kinetic 
Eqs. (fT2|) and p3|) . which is straightforward but cumbersome is presented in the Appendix. 

IV. SPIN ACCUMULATION 

To begin with the solution of the kinetic equation, which is derived in the Appendix, is 
used for the calculation of the spin accumulation. From Eq. ()A8|) . we obtain for the vector 
of the field-induced components of the spin-density matrix 

/(s) = -^[K X E)yn\eu)-, r, (20) 

■'^ ' ms^ ^ V2sr + 2a;2(2sr + l)' ^ ' 
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Fig. 1. Frequency dependence of the normalized 

non- vanishing component of the spin density 
matrix /^(w) for LOkpT — 0.1, 1, 2, 3, and 4 with 
jy{Q) = eETK/inh). 



Fig. 2. Time dependence of the normahzed 
non- vanishing component of the spin density 
matrix /^{t) for uikpT — 0.1, 5. The dashed hne 
is calculated from Eq. H23|l . 



with a = s + 1/t. According to Eq. (j2(J|) . the in-plane spin accumulation is calculated by 
a k integral over poles. Under the condition uj^t > 1, when one expects sharp resonances, 
the positions of these poles are given by 



1 

'2' 



[st) 



2,3 



-- ± 2iuJkT. 



(21) 



The resonance is most pronounced at zero temperature (T = 0) and when uj^t ^ 1. De- 
pending on the Rashba coupling constant a and on the carrier density, the resonance at 
2c(jfc^ {kp denotes the Fermi wave vector) is located in the THz regime. Switching from the 
Laplace to the frequency domain (s —iu) and considering zero temperature, we obtain 
from Eq. (jUJ 



K 



2iujt 



;2 ' 



(22) 



nhijjT [^ujI^ - [uj + i/rY] + 2iujl 
with ujkp = hKkF/m. It is a striking coincidence that the denominator in this expression for 
the spin accumulation completel y ag rees with the denominator in the spin-Hall conductivity 
calculated by Mishchenko et al. 12| • The associated spin excitation gives rise to resonances 
in the spin accumulation. An example is shown in Fig. 1 for some spin-orbit coupling 
parameters uJkpT. According to the fc-integral in Eq. (j2Uj) . the sharp peak is increasingly 
washed out with increasing temperature. The zero-frequency limit of the spin accumulation 
P{uj = 0) agrees with the result pubhshed by Edelstein |2(tI |. 

Applying an inverse Laplace transformation to Eq. (j^Uj) . the time evolution of the spin 
accumulation after the electric field is switched on can be determined. Numerical results are 
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shown in Fig. 2 (sohd hnes) and compared with the following analytical solution (dashed 
line) 

= ,Er- [l -exp [--) - 4;^exp (--jj , (23) 

valid under the condition uJ^pT ^ 1. As expected, for systems with a weak Rashba spin-orbit 
coupling, the steady state spin accumulation is only reached after a sufficiently long time. 
With decreasing oj^pT, weak oscillations are strongly suppressed. 

V. CHARGE CURRENT 

According to the consideration in Section III, the longitudinal current is calculated from 
the time derivative of the dipole operator [cf. Eq. (fT^ ]. Taking into account Eq. ()A12|) in 
the Appendix, we immediately obtain for the longitudinal charge current 

which can be used to calculate the frequency response at zero temperature (n' denotes the 
derivative dn{eu) / dsk). The result in the frequency domain 

. 2epT e'E 2zu^ 

Txh^ T^m ujT [AujI^ -{u: + i/rf] + 2iu:l^ ' ^ ' 

is composed of two contributions. The first one is expressed by the well-known Drude 
conductivity. The second one is due to the spin-orbit interaction and exhibits the same 
resonant denominator as the spin accumulation in Eq. (j22j) . At c<j = 0, Eq. (|25|) reproduces 
the result published in Ref. |8|- The measurement of the resonant longitudinal charge current 
contribution allows the determination of the Rashba coupling constant a. The alternative 
formulation in Eq. ()18p exactly reproduces Eq. ()24|) . 

The dynamical charge-Hall current in the Rashba split system has been studied recently 
I27I. The appearance of this effect is related to the k dependence of the scattering ma- 
trix elements. As this dependence is not taken into account in our approach, we obtain a 
vanishing charge Hall conductivity. 

VI. SPIN-HALL EFFECT 

To introduce the spin-Hall effect, let us first shortly recapitulate the main findings derived 
in the literature based on the "conventional" definition of the spin current in Eq. (fT^. 
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The treatment of the spin current within this framework reveals an anomalousness, which 
in our opinion was not duly noticed by many researchers. The approach predicts a non- 
vanishing stationary spin-Hall current that is independent of the electric field. Indeed, 
inserting Eq. ()A5|) into the expression in Eq. (jl9j) for the spin current, we obtain 

H K 

^.^(^) = ^E^'=^'(^-)' (26) 

k 

which leads to the constant x component of the spin-Hall current jy{uj) = —K£p/{2Tifi) in 
the frequency domain at zero temperature. The absence of any resonance indicates that 
there is no relationship between this fictitious spin-Hall current and the spin accumulation 
in Eq. (j^ . 

The field-induced spin-Hall current is calculated from its definition in Eq. (|T9|l and by 
taking into account Eq. ()A8|) . For zero temperature, we obtain for the frequency-dependent 
spin-Hall current 



which was recently derived by applying the Keldysh approach |12|. Here, the same res- 
onant denominator appears as in the longitudinal charge current [Eq. (|25|)] and the spin 
accumulation [eq. (EH)]- For a free electron gas (r oo), Eq. (j77j) simplifies to 



,2 



which was previously obtained by Erlingsson et al. Finally the steady state spin-Hall 

current {uj —>■ 0) of the clean Rashba model is given by the universal value 

m 

Although many authors confirmed these results, there remain some reservations. First of 
all, the definition of the "conventional" spin-Hall current led to considerable confusion and 
to serious doubts on its experimental relevance 22, 23 1. The main difficulty results from 
the fact that the spin is not a conserved quantity. As it has been claimed recently, a proper 
definition of the spin current requires a careful analysis of the torque density. It is assumed 
that this quantity may complement the above fictitious current to a conserved spin current. 
The prerequisite for such a construction is given, when the averaged spin-torque density 
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vanishes in the bulk |2^. It has been argued that this condition is fortunately satisfied for 
many spin models treated in the literature. This conserved spin current has the advantage 
that it can be measured via the spin accumulation to which it is related by the continuity 
equation. It has been pointed out that this quantity is straightforwardly calculated from 
the time derivative of the spin displacement (Eq. (5) in Ref. 22]). A firm foundation of 
this approach provides the definition of the spin current in Eq. ()17|) and the kinetic equation 
treated in Section II and solved in the Appendix. 

Based on the physically motivated definition of the spin current in Eq. ()17|) and using 
the analytical solutions of the kinetic equations presented in the Appendix, we shall restart 
the study of the spin-Hall current of the Rashba model. First, it is noted that we also get a 
spin-Hall current contribution that is independent of the electric field. From Eq. (jl7j) and 
()A11|) . we obtain a result 

.;W--f Z-M ...,^:!,.,^!) (30) 

that differs from Eq. ()26|) in many respects. First of all, this current contribution is closely 
related to the spin accumulation in Eq. ((201) • ^ ^^st evidence for this conclusion is the 
appearance of the same resonant denominator. However, this relationship goes even deeper 
and has a firm physical foundation, which becomes obvious by comparing the analytical 
solution for /o£ [Eq. ()A7|) ] with its counterpart for f^Q [Eq. ()A10|) ]. One solution is obtained 
from the other one by the replacement k, —ieEexde, where the derivative with respect 
to the energy applies to the charge density n. Transforming back this replacement to the 
spatial dependence, we obtain V ^ V + eEde [12], which gives the general recipe of a 
quasi-chemical potential to translate spatial inhomogeneities to internal field fiuctuations 
and vice versa. It is this connection and not the conservation of spins that establishes 
the close relationship between the field-induced spin accumulation and a field-independent 
spin current, which contrary to Eq. (PU)) vanishes in the stationary regime. The spin-Hall 
current in Eq. (jHUI) is interpreted as the response to the initial time evolution of the spin 
accumulation after the electric field is switched on. When the spin accumulation reaches its 
stationary value (cf. Fig. 2), the related spin-Hall current component disappears. 

Starting from the physically motivated definition of the spin current in Eq. (fTTj) and using 
the analytic solution in Eq. ()A13jl derived in the Appendix, we obtain the following general 
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result for the Laplace transformed spin-Hall current 

eEr fhKV ^ AtoKl + 2st) + 2a^ujl{l + 3st) - a^sr 

m J ^-^ 



h 



(31) 



^ [a^ST + 2cufc(2sr + 1)]' [a^sr + AuKst + 1)] ' 

which differs from all previous results in many respects. First, the denominator is composed 
of three factors, the zeros of which characterize spin eigenmodes. The in-plane spin precision 
is characterized by poles resulting from 1/ [a'^sr + AluK-st + 1)]^. Resonances of this kind 
appear both in the spin accumulation and in the charge current. In addition, there is the 
factor cr^sr + 4co'|(sr + 1) in the denominator of Eq. ()3ip that is associated with out-of-plane 
spin eigenmodes. The most striking discrepancy results from the k sum over the entire 
spin-orbit coupled Fermi sea, whereas only contributions from the Fermi surface are needed 
in charge transport problems. All states contribute to the time-dependent spin-Hall current. 
This astonishing result refers both to the spin current in [Eq. (IHUj)] and in [Eq. (jS])]. 

To clarify the origin of this peculiarity, we note 
that the spin current is not due to displacements 
of carriers but induced by the change of the mag- 
netic moment. This situation is quite similar to 
the diamagnetism in normal metals, which is also 
determined by all states in the entire Brillouin zone 




28|. Finally, we point out that the k integral in 



Eq. (j3Tj) leads to a logarithmic singularity, the due 
treatment of which requires a careful consideration 

Fig. 3. Time dependence of the normalized kinetic equations in the limit UJf, 0. 

spm-Hall current for ujkpT = 5 (thick sohd Applying an integration by parts, we obtain an 
hne), 1 (thick dashed hne), and 0.1 (thin soUd 

hne) at T = 0. equivalent expression for the spin-Hall current, in 

which the logarithmic contributions are singled out and in which at zero temperature the 
k integral is artificially fixed at the Fermi surface: 



eEr 
2h m 



2ujU2st + 1] 



[a^sT + 2ul{2sT + 1)] 



(32) 



:i + sr)(8sr + 3) 
2(2sr + 1) 



In 



1 + 



2u;l{2sT + 1) 



{sry + 15sr + 6 
4(sr + 1) 



In 



1 + 



From this equation, it is concluded that in the steady state (s — » 0) the spin-Hall current 
vanishes for the impure Rashba model of two-dimensional electrons. This finding agrees 
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completely with former conclusions derived from the "conventional" definition of the spin- 
Hall current 

Applying an inverse Laplace transformation, the time evolution of jy can be studied. 
Fig. 3 shows an example for the time-dependence of the spin-Hall current, which is induced 
by switching on the electric field at t = 0. Depending on the coupling parameter uj^pT, 
strong oscillations of the spin-Hall current initially develop, which are completely damped 
out after a couple of scattering times. 

If time-dependent electric fields are applied, the spin-Hall conductivity cr^/f becomes 
nonzero. The frequency-dependent spin-Hall conductivity is obtained from Eq. (jH^ by an 
analytic continuation {s — > —iuo). Numerical results for the real and imaginary part of 
are shown in Fig. 4 and 5 by thick solid lines. We focus on the zero-temperature case and 
compare with previous results obtained from the "conventional" definition of the spin-Hall 



current [dashed lines as calculated from Eq. (jlTj)]. Both approaches predict a sharp 
resonance in the ac spin-Hall conductivity, when the condition UkpT ^ 1 is satisfied. The 
enhancement of the spin-Hall current appears at = 2ujkp- As seen from Fig. 4, remnants 
of this feature survive even in the case oJkpT < 1. In the limit uikpT ^ 1, both approaches 
agree and result in 



eE 



2'Kfi V 1 — iujT 



(33) 



,x 10 



Re CT,f//cr(i 
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Fig. 4. Semilogarithmic plot of the real and 
imaginary parts of the normalized spin-Hall 
conductivity as a function of ujt for UkpT — 0.1 
(thick solid line). The dotted (dash-dotted) line 
is calculated from Eq. [Eq. ^^]. o-q is 
given by e/(2nh). 




Fig. 5. Real and imaginary parts of the 
normalized spin-Hall conductivity as a function 
of WT for u!kpT — 1 (thick solid line). The 
dashed line is calculated from Eq. (|27|l . 
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which is plotted by the dash-dotted hne in Fig. 4. From Eq. ()33p . the non-analytic behavior 
of the spin-Hall current becomes obvious because the approximation fails in predicting zero 
spin-Hall current in the limit {u 0), which requires the general result in Eq. (|31|). 

Comparing the solid and dashed lines in Figs. 4 and 5, it is tempting to conclude that there 
is no qualitative difference between the results of both approaches. That this conclusion is 
only partly true shows a consideration of clean samples (r oo). In this case, our approach 
becomes completely exact and we obtain 

= 8^(2 " ^ (4<-c.2). + 2(4<-.2) - 2I' 

which gives a logarithmic divergency at vanishing frequency 

Most previous approaches predict in this case a finite universal spin-Hall conductivity (Jsh = 
In contrast, based on a physically motivated definition of the spin-Hall 
current and on an exact procedure, we obtain neither zero nor an universal value but a 
logarithmic divergency for the spin-Hall conductivity at zero frequency. 



Unfortunately, recent calculations [22| of the physical spin-Hall conductivity in the clean 
limit of the non-interacting Rashba model do not agree with our exact result in Eq. (jH^jl . 
The authors started from the same physically motivated definition of the spin-Hall current 
(Eq. (5) in Ref. and calculated jy via the sourceless continuity equation. However, 

strictly speaking, the spin current, which is directly calculated from the k derivative of the 
averaged distribution function f{k,K, \ t), is not related to source but to vortex fields. To 
illustrate the situation, let us treat the set of kinetic equations for E = and to first- 
order in K. These equations contain two curl contributions namely ~ [K x K]f{k, k, \ s) and 
~ [kx f{k, K, \ s)], which do not enter any continuity equation. The first vector describes the 
coupling between charge and spin degrees of freedom and gives rise to the field-independent 
spin-Hall current. If an electric field is switched on, the replacement k — > k — ieEde 
generates a contribution ~ —ie[K x E]def{k, n \ s), which is responsible for the appearance 
of the field-induced spin accumulation. For the mean values f{e \ s) = f{k, k | s) and 
j{e I s) = isd^ ® f{k^ n \ s) /2, which are used to express the spin accumulation and 
the field- independent spin-Hall current, respectively, we obtain a relationship of the form 

'-f\e I s)n = eE,j^{e \ s)n\ (36) 
14 



which is confirmed by our approach [compare Eq. ()2U|) and Eq. ()30|) ]. This equation is not 
vahd for the "conventional" field- independent spin-Hall current in Eq. The interrelation 
in Eq. (j36|) is based on the assumption that the field effects can be accounted for by a quasi- 
chemical potential in the electronic part of the density matrix. For a = 0, it is known 
that this supposition leads to the Einstein relation between the mobility and the diffusion 
coefficient. Recently, this hypothesis has also been accepted for systems with spin-orbit 
interaction |l2, 29, 3y|, although its application becomes more subtle due to the spin splitting 
of the energy bands. An analogous use of quasi-chemical potentials for the treatment of field 
effects on the spin density requires further justification. 

The appearance of the logarithmic dependence of the ac spin-Hall conductivity in the 
disordered two-dimensional Rashba model remind us very much of the logarithmic quantum 



corrections in the theory of weak localization 



3l|. Based on scaling arguments, it was shown 



that the ac conductivity of a disordered two-dimensional electron gas exhibits a logarithmic 
divergency in the zero-frequency limit [cr ~ In(co'r)]. Our result for the spin-Hall conductivity 
in Eq. (jHKjl is comparable to this dependence observed in the completely other field of weak 
localization. 



VII. SUMMARY 

Based on the kinetic equations for the spin-density matrix of the two-dimensional Rashba 
model, we treated the electric-field-induced spin accumulation and spin transport in the 
linear response regime. At zero temperature, the frequency dependence of both the spin 
accumulation and the longitudinal charge- carrier transport exhibit a sharp resonance at 
u = 2ujkp {ujkp = akp/h, with a being the spin-orbit coupling constant), which is due to 
eigenmodes of spin excitations. The measurement of this resonance should be possible. It 
allows the determination of the Rashba coupling constant a. Similar resonances are expected 
to appear also in the fc-cubed Rashba model for 2D holes and the Luttinger model for 3D 
holes. 

For the charge-carrier transport, there are two completely equivalent procedures to calcu- 
late the conductance. The same results are obtained by starting either from the symmetrized 
product of the density and velocity operators or the time derivative of the dipole moment. 
Unfortunately, this equivalence no longer holds for the spin transport. The spin is not a 
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conserved quantity. Therefore, it is a serious problem to chose a proper definition for the 
spin current that does not lose its physical foundation. Most researchers preferred a defini- 
tion of the spin-Hall conductivity, which revealed rather unconventional properties so that 
serious doubts arose on its experimental relevance. Recently, a physically motivated defi- 
nition of the spin current has been suggested [22] that resolved a number of difficulties of 
former approaches. It has been argued that the proper effective spin current is inevitably 
defined as the time derivative of the spin displacement. Applying this definition, results for 
the spin-Hall conductivity are obtained that are drastically different from previous findings. 
First, the approach predicts a field-independent spin-Hall current that refiects the initial 
variation of the spin accumulation after the electric field is switched on. Contrary to previ- 
ous results, this specific spin current contribution disappears in the steady state. Its physical 
origin is due to the initial time evolution of the spin polarization. Furthermore, for a clean 
two-dimensional Rashba model, we obtain a spin-Hall conductivity that exhibits a logarith- 
mic dependence at low-frequencies and not an universal constant value. This observation, 
which is an exact result, reminds us on the well-known frequency-dependent conductivity of 
a disordered two-dimensional system in the theory of weak localization. 
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APPENDIX A: SOLUTION OF THE KINETIC EQUATIONS 

The kinetic Eqs. fll2|l and (fT!?|l are solved by an perturbational approach with respect to 
E and k. This calculation exploits the formal exact solutions of these equations given by 

o-r - 2a;fc X r -F4a;fc(c5fc ■ r)/cr 

^ - ^^T^, ' ^^^^ 

with a = s + 1/t, r = R + f /t and 

— * 

ih , , ^ ih , , eE 

R=-{k- k)f - -{K X K)f - —Vkf 

mm h 

+-—fhuJk jz—f- A2 
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For the angle-averaged spin-density matrix, we obtain 



for their x, y components, while for the z component it follows 



-J- aR-2u}ky<R .... 
P = . A4 



The lowest-order solutions in £^ = and k = (the corresponding elements of the density 
matrix are denoted by /oo and /oo) are easiliy obtained 

/oo = ^ ^ = /oo; /oo = —fi^k — , /oo = 0, (A5) 
s s 

where n' denotes the derivative with respect to Sk- Both quantities /oo and /oo do not depend 
on time and are therefore conserved under the condition of thermodynamic equilibrium. For 
the derivation of this result it was necessary to consider the spin-orbit coupling in the 
collision integral [the second and third term on the right-hand side of Eq. (jl3|) ]. Next, let us 
calculate the lowest-order correction due to the electric field [E 7^ 0, k = 0). Taking into 
account Eqs. (fT^ and (fT^ together with Eq. ()A3|) . we obtain 

foE=- — —n', 7^=0, (A6) 
as m 

foE = uJkieEk^) + n r — — — —- , A7 

mas [a^sT + 2uJk[2sr + 1)\ 

I- <^E i , ,N, 2aru;ln' 1 , . 

as I a^sT + 2u!l{2sT + 1) J 

where we used the abbreviation uje = eh{K x E)/m. Next, we calculate /^o and f^a by 

collecting the first-order contributions in n and by setting E = Q. Neglecting corrections, 

which are of the order h^K/msy., we obtain 

ih 

fno = (K-fc)n, (A9) 

mas 

F' /, x^^^' . o-o)^ - 2u3fc X o;^ + 4(a3fc ■ a;^)a3fc/c^ ..,„^ 
/«o = -«t^fc fe ■ K mr ^ — — —- AlO 

mas [a^sT + 2uJk{2sT + 1)J 

1 y 2aTujln \ 

as a'^ST + 2u!l[2sT + 1) J 
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There is an interesting symmetry between the vectors /^o and foE- From Eq. ()A10|1 . 
the electric-field-induced contribution in Eq. ()A7jl is obtained by the replacement k, — *■ 
—ieEexds, where the derivative refers specifically to the charge density n. Finally, we need 
the first-order corrections in k, and E of the angle-averaged component of the density matrix, 
which are expressed by 

f.E = -\—PoE + [2{eny-n] , A12 

s [ m mas J 
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